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We demonstrate the onset of strong on-site localization in a one- dimensional many-particle system. 
The localization is obtained by constructing, in an explicit form, a bounded sequence of on-site 
energies that eliminates resonant hopping between both nearest and remote sites. This sequence 
leads to quasi-exponential decay of the single-particle transition amplitude. It also leads to strong 
localization of stationary many-particle states in a finite-length chain. For an infinite chain, we 
instead study the time during which all many-particle states remain strongly localized. We show 
that, for any number of particles, this time exceeds the reciprocal frequency of nearest-neighbor 
hopping by a factor ~ 10 5 already for a moderate bandwidth of on-site energies. The proposed 
energy sequence is robust with respect to small errors. The formulation applies to fermions as 
well as perpetually coupled qubits. The results show viability of quantum computing with time- 
independent qubit coupling. 

PACS numbers: 03.67.Lx,72.15.Rn,75.10.Pq,73.23.-b 



I. INTRODUCTION 

Disorder-induced localization has been one of the cen- 
tral problems of condensed matter physics, starting with 
the Anderson paper Tj. It has also attracted much at- 
tention in otherphysical contexts, quantum chaos being 
an example 0, . Recent interest in quantum comput- 
ing has further emphasized its importance and allowed 
looking at it from a somewhat different perspective. 

In many proposed physical implementations of a quan- 
tum computer (QC) the qubit-qubit interaction is not 
turned off 0, 0, 0, S 13 H3 • The interaction may lead 
to excitation hopping from one qubit to another. How- 
ever, control and measurement should be presumably 
performed on individual qubits. Therefore it is essential 
to prevent excitation transfer between operations. This 
makes localization a prerequisite for quantum comput- 
ing with perpetually coupled qubits. Several approaches 
to quantum computing with perpetually coupled qubits 
have been proposed recently [ll], ll3> D3 

In a multi-excitation system like a system of interact- 
ing electrons or a QC, inter-site (or inter-qubit) transi- 
tions are a many-body effect, they involve several exci- 
tations. In spite of the broad interest in the problem 
of many-particle localization, only a limited number of 
analytical results has been obtained 0]. Numerical re- 
sults are also limited: classical computers do not allow 
studying a large number of particles, because the Hilbcrt 
space is exponentially large. On the other hand, QC's 
with perpetually coupled qubits provide a unique means 
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for investigating localization in a controllable setting. 

In this paper we study strong on-site many-particle lo- 
calization. It implies that each particle (or excitation) is 
nearly completely confined to one site. This is a stronger 
condition than just exponential decay of the wave func- 
tion, and it is this condition that must be met in a QC. 

A well-known argument suggests that it is hard 
to strongly localize a disordered many-particle system 
where the on-site energies are random and uniformly dis- 
tributed within a finite- width band 0] . Indeed, consider 
a state where particles occupy N sites. For short-range 
hopping, it is directly coupled to oc N other TV-particle 
states. With probability oc N one of them will be in 
resonance with the initial state. For large iV this leads 
to state hybridization over time ~ J -1 , where J is the 
inter-site hopping integral (we set H = I). 

In a QC, the quantity J is determined by the qubit- 
qubit interaction and usually characterizes the rate of 
two-qubit operations. At the same time, the qubit ener- 
gies are not random and often can be individually con- 
trolled. It is this control that makes QC's advantageous 
for studying many-particle localization, as it becomes 
possible to construct a "disordered" energy sequence site 
by site. 

Localization can be considered from two points of view. 
One is based on the analysis of stationary states of a 
many-particle system. The other is based on studying 
the system dynamics. Consider a state with particles 
occupying a given set of sites, which is called an on-site 
state (or a quantum register) . As a result of hopping it 
can hybridize with another on-site state with nearly the 
same energy. We will study the time it takes for resonant 
hybridization to happen, which we call the localization 
lifetime ii oc . 
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In a QC all states have a finite coherence time due 
to coupling to the environment and external noise. For 
successful QC operation, derealization should not occur 
during this time. For most of the proposed models of a 
QC, the coherence time is < 10 5 J -1 . Therefore it is suf- 
ficient to have the localization lifetime > 10 5 J~ 1 . Such 
lifetime-based formulation of the many-particle localiza- 
tion problem is relevant to condensed-matter systems as 
well, because of finite decay and decoherence times of 
quasiparticles for nonzero temperatures. 

Here we construct a bounded sequence of on-site en- 
ergies in a ID chain and show that it leads to a long 
localization lifetime. We provide evidence that it also 
leads to strong localization of many-particle stationary 
states in sections of the chain with length up to 12 sites. 
In a sense, this is an explicit construction of an efficiently 
localizing on-site disorder. 

In a QC, on-site excitation energies are interlevel dis- 
tances of the qubits. They can often be individually con- 
trolled, which makes it possible to construct an arbitrary 
energy sequence. However, since the qubit tuning range 
is limited, so should be the energy bandwidth. A smaller 
bandwidth leads also to a higher speed of quantum gate 
operations, particularly if they involve changing qubit en- 
ergies 16]. Of course, condensed-matter systems always 
have bounded bandwidth as well. 

In condensed-matter physics, localization by "con- 
trolled" disorder has been studied in depth and many 
interesting results have been obtained in the context of 
incommensurate periodic potentials, see 

Refs.mEm 

and papers cited therein. In contrast to this work, 
we are interested in strong on-site localization, and not 
only for single-, but in the first place, for many-particle 
states. The "potential" that we propose is not quasi- 
periodic, and our analytical techniques, including the 
time-dependent formulation for many-body systems, are 
different from the methods developed for quasi-periodic 
potentials. 

To strongly localize one particle, the difference between 
excitation energies on neighboring sites should be much 
larger than the hopping integral J. However, even for 
nearest neighbor coupling, the energies of remote sites 
should also differ to prevent transitions between them 
via intermediate (nonresonant) sites. The further away 
the sites are, the smaller their energy difference can be. 
This idea is implemented in our energy sequence. As a 
consequence, the single-particle transition amplitude dis- 
plays nearly exponential decay with distance. We show 
that the decay exponent weakly depends on site and find 
rigorous bounds on its value. 

For many-particle localization one has to suppress not 
only single-particle, but also combined resonances, which 
involve simultaneous transitions of several interacting ex- 
citations. There is no known way to eliminate all such 
resonances in an infinite system. However, with the in- 
creasing "order" of the transition, i.e., the number of 
involved excitations and/or intermediate virtual states, 
the effective hopping integral may quickly fall off, lead- 



ing to an increase of the transition time. Then to obtain 
a desired lifetime of a localized state it is sufficient to 
eliminate resonances up to a certain order. We explicitly 
show how to do it up to fifth order, for our sequence. 

In a real system it will be possible to tune the energies 
only with certain precision. We study how errors in the 
energies affect localization and show that our sequence 
is stable with respect to small errors. We also demon- 
strate that, in terms of strong on-site localization, even 
for a small chain the constructed energy sequence is much 
better than fully random on-site energies with the same 
overall bandwidth. 

The paper is organized as follows. In Sec. II we discuss 
the Hamiltonian of coupled qubits and introduce a phys- 
ically motivated one-parameter sequence of on-site en- 
ergies. In Sec. Ill one-particle localization of stationary 
states is studied and quasi-exponential decay of the tran- 
sition amplitude is demonstrated. A rigorous proof of 
such decay is provided in Appendix A. The relevant scal- 
ing properties are analyzed in Appendix B. In Sec. IV the 
inverse participation ratio is calculated for many-particle 
states in a section of a ID chain. This ratio can be made 
very close to one in a broad range of the parameter of the 
on-site energy sequence, but it also displays sharp reso- 
nant peaks as a function of this parameter. These peaks 
are discussed in Appendix C. The lifetime of localized 
states is considered in Sec. V. A minor modification of 
the energy sequence allows one to open a gap in the spec- 
trum of combined many-excitation transitions up to 5th 
order, which is sufficient for extremely long localization 
lifetime. The role of errors in the on-site energies is stud- 
ied, and robustness of the results with respect to these 
errors is demonstrated in Sec. VI. In Sec. VII a highly 
symmetric period-doubling sequence of on-site energies 
is analyzed along with a sequence of uncorrelated on-site 
energies. Both are by far inferior, in terms of localiza- 
tion, to the sequence discussed in Sees. II-V. Sec. VIII 
contains concluding remarks. 



II. THE MODEL 

The problem of localization can be formulated in a 
similar way for one-dimensional systems of interacting 
fcrmions and spins. The formulation also applies to 
qubits during the time when no gate operations and 
measurements are performed, i.e., there are no time- 
dependent fields that would modulate the qubits. 

The relation between spin and qubit systems is sim- 
ple: qubits are two-level systems, and therefore can be 
described by S = 1/2 spins in a magnetic field. Then the 
excitation energy of a qubit becomes the Zeeman energy 
of a spin, whereas the qubit-qubit interaction becomes 
the exchange spin coupling. Note that the physical in- 
teraction itself may be of a totally different nature, e.g., 
electric dipolar or quadrupolar. 

For many proposed realizations of QC's the qubit ex- 
citation energies are large compared to the interaction. 
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Then the qubit-qubit interaction is described by the spin- 
coupling Hamiltonian of the form 
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Here, n, m enumerate sites in the ID spin chain, z is 
the direction of the effective magnetic field, and are 
the interaction parameters (/i = x, y, z). In mapping the 
qubit interaction on the exchange coupling we kept only 
those terms which, in the Heisenberg representation, do 
not oscillate at qubit transition frequencies. This is why 
we have set = J^m- ^ n the case °f qubits, the terms 
S£S* = (1/2)(S+S-+S-S+) lead to excitation 
transfer between the qubits n, m provided their energies 
are close. We note that the spin interaction of the form 
Q conserves the number of excitations in the system. 

It is convenient to map the spin system onto a system 
of spinless fermions via the Jordan- Wigner transforma- 
tion [21|. For nearest neighbor coupling, the Hamiltonian 
of the fermion system becomes 



H — Hn + Hi 



(2) 



Ho = e n a} n a n + ^ J ^^(<4a n+ i + a\ l+l a n ), 

* — 

Here, aj^dn are the fermion creation and annihilation 
operators; the presence of a fermion on site n corre- 
sponds to the nth spin being excited. The parameter 
J = J xx n+l is the fermion hopping integral. The parame- 
ter J A = J„ z n+ i gives the interaction energy of fermions 
on neighboring sites. If the coupling of the underlying 
spins is isotropic, we have A = 1. The on-site fermion 
energies e n are the Zeeman energies of the spins (exci- 
tation energies of the qubits) counted off from the char- 
acteristic central energy which is the same for all spins 
inside the chain. For concreteness we set J, A > 0. 

Localization, and in particular weak localization, is of- 
ten described in terms of the decay of the wave functions 
of stationary states at large distances. In contrast, here 
we are interested in strong on-site localization. It is de- 
termined by the short-range behavior of stationary states 
and corresponds to confinement of each particle to one 
or maybe a few neighboring sites. 

Strong localization can be conveniently characterized 
by the inverse participation ratio (IPR) , which shows over 
how many sites the wave function spreads. For an N- 
particle wave function \iI>n\) (A enumerates the station- 
ary states) the IPR is given by the expression 
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where |0) is the vacuum state, 
we will sometimes use the notation 

J „t 



a ni \ipN\)\ 

(3) 

In what follows 
|$(fci,fc 2 ,...)) - 



0) for the on-site wave function (quantum reg- 
ister) in which sites k\ , ki , . . . are occupied and other sites 
are empty. The quantity Q is also sometimes called the 



number of participating components, it shows how many 
matrix elements ($(fci, fe, ■ ■ OIV'jva) are substantial, i.e., 

of order of their maximal value with respect to ki, fe, 

For fully localized states Inx = 1- Strong on-site lo- 
calization means that Inx is close to 1 for all states A. 
In this case both the average IPR 



NX 

and In max = m&xxlNX are close to one. Here, Cnx is 
the total number of TV-particle states; for an L-site chain 
C N x = L\/N\{L - N)\. Smallness of (I N ) - 1 is a weaker 
condition, it is an indication of strong localization of most 
of the states. 

In the opposite limit of extended states we have 
{Inx) ~ Cnx 3> 1. A simple example is the case of one 
particle in an open chain (a chain with free boundaries) 
with e„ =const. The wave functions of the particle are 
sinusoidal, and for an L-site chain (7i) = 2(L + l)/3. The 
mean IPR sharply increases with the number of particles 
N, for N < L/2. 



A. The on-site energy sequence 

Localization requires that the on-site energies e n be 
tuned away from each other. The strategy for choosing 
the sequence of s n while keeping the overall bandwidth 
of the energy spectrum finite can be as follows. First, 
we separate the energies of nearest neighbors by splitting 
e n into two subbands, with even and odd n, respectively. 
The distance h between the subbands should significantly 
exceed the hopping integral J. Then each subband is 
further split into two subbands in order to detune next 
nearest neighbors. The splitting between these subbands 
can be less than h. This is because next-nearest-neighbor 
hopping occurs via virtual transitions to a nonresonant 
nearest-neighbor site, and therefore the effective hopping 
integral is ~ J 2 /h. The procedure of band splitting 
should be continued, and higher-order splitting can be 
smaller and smaller. 

We now introduce a simple sequence of s n that im- 
plements the structure described above. Except for the 
energy scaling factor h, this sequence is characterized by 
one dimensionless parameter a. As we show, it can al- 
ready be efficient in terms of strong localization. For a 
semi-infinite chain with n > 1 we set 
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[n/k\ a k-l 



(4) 



([•J is the integer part). 

Sequence Q does not have any simple symmetry. For 
example, it is not self-similar: the subband widths do 
not scale with the distance between the sites that be- 
long to the same subband (cf. Sec. VII). Nor is sequence 
Q quasi-periodic. However, the coefficients at any given 
power a q are repeated with period 2(q + 1). This impor- 
tant property is essential for obtaining analytical results, 
see Appendix A. 
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The energy spectrum Q is illustrated in Fig. ^ The 
left panel gives the energies of the first 50 sites. It is 
seen that the sites with close e n are spatially separated, 
whereas the sites with close n are separated energeti- 
cally. The multisubband structure of the spectrum is 
clearly seen in the right panel. For small a, the two ma- 
jor subbands have width w ah and are separated by ~ h. 
The splitting of higher-order subbands is proportional to 
higher powers of a. As a increases the subbands start 
overlapping, and for a > 0.4 separation between the sub- 
bands disappears. 

As a result of the low symmetry, different subbands in 
the right panel of Fig^have different numbers of points, 
i.e., e n are not evenly distributed among the subbands. 
This turns out to be important for strong many-particle 
localization. The case of a symmetric sequence is dis- 
cussed in Sec. VII. 
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FIG. 1: (color online) The energies e n /h for a — 0.3. The 
left panel shows e n /h for the sites n = 1,2, ...,50. Sites 
with close on-site energies are spatially separated. The right 
panel shows e„/h for a much longer array, n = 1, . . . , 2000. 
The energy spectrum displays a multisubband structure, with 
clearly identifiable 16 subbands in this case. 

An important advantageous feature of sequence for 
quantum computing is that it is convenient for perform- 
ing gate operations. For single-qubit gates, a single par- 
ticular radiation frequency can be used to resonantly ex- 
cite different qubits. It has to be chosen near the average 
single-qubit transition frequency (which corresponds to 
e = 0). Then qubits can be selectively excited by tun- 
ing them to this frequency. The transition frequency of 
the qubit depends on whether neighboring qubits are ex- 
cited. This can be used for implementing a CNOT gate. 
Alternatively, neighboring qubits can be tuned in reso- 
nance with each other, which will lead to a two-qubit 
excitation swap |l6j . 

III. SINGLE-PARTICLE LOCALIZATION: 
STATIONARY STATES 

A. The transition amplitude 

In a ID system with random on-site energies all single- 
particle stationary states are localized, even for weak dis- 
order, and exponentially decay at large distances. Al- 



though sequence is not random, the transition ampli- 
tudes also display quasi-exponential decay at large dis- 
tances provided J <C ah, as follows from the results of 
Appendix A. In this paper we are primarily interested 
in the short-range behavior. It turns out that a parti- 
cle is confined much stronger in the case of sequence Q 
than in the case of random on-site energies distributed 
within the same energy band, see Sec. VII. The confine- 
ment quickly strengthens with the increasing parameter 
a once a exceeds a certain threshold value a t h- 

Spatial decay of single-particle stationary states can 
be characterized by the amplitude of a particle transition 
from site n to site n + m. To the lowest order in J/h it 
has the form 

nm 
J/2(e n -s n+k . (5) 

For sequence (@J in the limit of small a the energy 
difference \e n + m — £ n \ can be approximated by its leading 
term, so it is ~ h for odd m and ~ ah for odd m/2. In 
general, the larger is m the higher may be the order in a 
of the leading term in |e„+ m — e n \. 

The asymptotic behavior of the function K n (m) for 
small a and large m can be studied rigorously. The anal- 
ysis is based on some results of number theory. It is 
given in Appendix A. It shows that K n (m) decays with 
the distance m quasi-exponcntially, 

K n (m)=a- u W(J/2h)W. (6) 

The exponent v is determined by limlog-ftr„(TO)/m for 
m — ► oc and depends on n . The values of v are bound to 
a narrow region centered at v = 1, with 0.89 < v < 1.19. 
For estimates one can use v = 1, i.e., 

K n {m)7aK\ m \ K = J/2ah. 

The decay length of the transition amplitude is then 
X/\\nK\. 

The numerical values of v for different n and m are 
shown in Fig. [5] They were obtained by keeping the 
leading term with respect to a in each energy difference 
£n — £n+k with 1 < k < m. The data are in excellent 
agreement with the asymptotic theory. 

We have also studied decay in the "opposite" direc- 
tion, i.e., for negative m in Eq. (JSJ) (however, m + n > 1). 
The asymptotic expression © applies in this case, pro- 
vided n, |m| 3> 1. The numerical data also agree with 
the theory. 

Equation © gives the tail of the transition amplitude 
for J/2ha -c 1. It does not immediately describe strong 
single-particle localization, which is determined by the 
short-range behavior of the wave function. However, one 
may expect that strong localization should occur when a 
becomes much larger than a typical threshold value ath = 
J/2h. In fact, Eq. @) describes the transition amplitude 
only when ath C a <C 1. We note that the inequality 
a t h < 0.4, which is necessary to avoid overlapping of 
the subbands, is satisfied already when the ratio of the 
energy bandwidth to the hopping integral hj J exceeds 
the comparatively small value 1.3. 



5 




<>1> 



FIG. 2: (color online). The exponent v of the a-dependence 
of the transition amplitude K n (m) for the efficient distance 
m — 200 (crosses) and m = 1000 (circles) as a function of site 
number n. The dashed lines show the analytical limits on v. 



B. The inverse participation ratio 

A quantitative indication of strong localization of 
single-particle stationary states is that 1\\ -1 < 1 for all 
states A. Numerical results on (ii) and I\ max as functions 
of the energy spectrum parameter a for two values of the 
scaled bandwidth h/J are shown in Fig. EI They were ob- 
tained by diagonalizing the Hamiltonian @ numerically. 
The data refer to open chains of three different lengths 
L, with the first site being always n = 1 in Eqs. J2J, I0J- 
The sum over n in the terms oc J, JA in Eq. @ ran from 
n = 1 to n = L — 1. 

In the fermion Hamiltonian J3J, the energies e n differ 
from the Zeeman energies of spins by — JA inside the 
chain. On the boundaries this shift is — JA/2. To make 
the Zeeman energies equal to e„+const, with the con- 
stant being the same for all spins, we added the term 
-| JA(Sl + SI) to the spin Hamiltonian JTJ. Then the 
numerical results for a finite open chain equally apply to 
both spin and fermion systems. 

In the limit a — > 0, the energies of single-particle 
stationary states form two bands centered at h/2 for 
even n and — h/2 for odd n. For h ^> J the widths 
of these bands are ~ J 2 /h. The system is equivalent 
to two weakly coupled translationally-symmetric open 
chains; the band wave functions are sinusoidal, which 
gives (ii) = (L + 2)/3. This agrees with the value of (7i) 
for a — > in Fig. 01 

For nonzero a the on-site level detuning (jJJ breaks 
translational symmetry. As a increases, the bands at 
±h/2 are split, and more and more subbands are resolved 
in the energy spectrum. Respectively, decreases. 
It sharply drops to 1 in a narrow region, which can 
be conditionally associated with a transition to strong 
localization. The center of the transition region gives 
the threshold value ath of the parameter a. It appears 
to be independent of the chain length L. The estimate 
ath = J/2h from Eq. © is in good agreement with the 
numerical data for different h/J. 

When ath -C a •< 1, all stationary states are strongly 
localized. Tails of the wave functions are small and lim- 
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FIG. 3: (color online). The mean single-particle inverse par- 
ticipation ratio (Ji) vs. a for h/J — 20 (upper panel) and 
h/J — 10 (lower panel). The data refer to three values of the 
chain length L. The vertical dashed lines show the analytical 
estimate for the threshold of strong localization. The insets 
show the maximal IPR over all eigenstates, I\ max = max* I\\ . 
It sharply decreases with the increasing a. The peak of Ji max 
for L = 300 near a — 0.1 is due to the boundary. Near 
the minimum over a, we have ii m ax ~ 1.01 for h/J — 10, 
and /imax ~ 1.003 for h/J = 20. This demonstrates strong 
single-particle localization. 



ited mostly to nearest neighbors, which leads to 



hx - 1 » J 2 /h 2 



(7) 



to lowest order in J/h,a. 

For a > 0.4 the IPR increases above its minimal value. 
This happens because the major bands of e n centered at 
±h/2 start overlapping. 

The minimum of the IPR over a is broad for large 
h/J. Near the minimum the numerical data in Fig-yjare 
in good agreement with the estimate (Q. The agreement 
becomes better with increasing h/J. 

The insets in Fig. [21 show that the IPR as a function 
of a can have narrow resonant peaks. In the presented 
data they occur for the chain of length L — 300. The 
peaks are seen only in Ii ma x, whereas (Ji) remains close 
to 1. This indicates that only a few on-site states are 
hybridized with each other. 

The underlying resonance results from a different 
hopping-induced shift of the energy levels at the chain 
edges compared to the bulk. The analysis of the wave 
functions shows that the peak corresponds to a resonance 
between sites 300 and 296. Because of the hopping, the 
energy of site 300 is shifted by w (J/2) 2 /h, whereas for 
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site 296 this shift is « J 2 /2h. The difference of e n for 
these sites is ~ ha 3 for a <C 1. Then the peak should oc- 
cur at a ~ (J/2/i) 2 / 3 , in good agreement with the data. 
The effective hopping integral between the two sites is 
determined by virtual transitions via intermediate sites, 
it is ~ J 4 /16ah 3 . Therefore the width of the peak with 
respect to a should be oc (J/h) 2 , also in agreement with 
the data. 

In Appendix B we outline another way of looking at 
the effect of the parameter a on localization. Specifically, 
we study the scaling relations between a and J/h that 
follow from the condition that the IPR takes on a given 
value. 



IV. MANY PARTICLE LOCALIZATION: 
STATIONARY STATES 

A. Many-particle hopping 

The many-particle localization problem is more com- 
plicated than the single-particle one. When the param- 
eter of the inter-particle interaction A ^ (i) the en- 
ergy levels e n are shifted depending on the occupation of 
neighboring states, e„ — > s n , (ii) there occur combina- 
tional many-particle resonances e ni + . . . + e nk w i mi + 
. . . + i mk , and as a result, (iii) there occur interaction- 
induced many-particle transitions that may be resonant 
even though single-particle resonances have been sup- 
pressed. Such transitions may lead to derealization. 

In contrast, the case A = corresponds to the XY- 
type coupling between the underlying spins. In this case 
the single-particle results apply to the many-particle sys- 
tem. 

To analyze many-particle effects, it is convenient to 
change from , a n to new creation and annihilation op- 
erators b n that diagonalize the single-particle part Hq 
of the Hamiltonian J5J, a n = J2k U n kbk- The unitary 
matrix U is the solution of the equation 



(^o)nm — ^n^nm 



e' s 



(8) 



Here, £' n are the exact single-particle energies, 



£/t,ff t/ = V :: 



' bH 

n u n u n- 



(9) 



For a 3> ath and J <C h, when single-particle states 
are strongly localized, the energies £' n are close to the 
on-site energies e n . To leading order in J/h, a we have 



2 r 



J_ 

2h 



(-!)« + -(-!) 



L«/2J , 



(10) 



The major term in the right-hand side corresponds sim- 
ply to renormalization of the characteristic bandwidth of 
on-site energies h — ► h + J 2 /2h. 



In terms of the operators b n , b^ the interaction part of 
the Hamiltonian is 

tfHiU = JAj2Vk lk2 k 3 kAA b k3 b k 4 , (11) 
where the sum runs over k\ 2 3,4, and 



Vfeifc 2 fc 3 fc4 — ^ U* kl U* +lk2 U p+ ik 3 U } 



pk 4 ■ 



(12) 



The Hamiltonian i|ll|) describes the interaction of the 
exact single-particle excitations. 

If the single-particle states are all strongly localized, 
the off-diagonal matrix elements U n k are small. They 
are determined by the decay of the wave functions, and 
therefore fall off exponentially with increasing \k — n\. 
From Eq. © we have U nk ~ ifl fe -™l for |fc - n\ » 1. At 
the same time, the diagonal matrix element is U nn ~ 1. 

Therefore, for strong single-particle localization the 
major terms in the matrix Vk 1 k 2 k 3 k i are those with x = 0, 
where 

x = min(|jfei-p| + |fe-i>-l| + |A!3-p-l| + |A!4-p|). (13) 
p 

These terms lead to an energy shift cx J A of the states de- 
pending on the number of particles on neighboring sites. 

The meaning of the parameter xr (|13|) can be under- 
stood by noticing that the terms oc V kl k 2 k 3 k 4 m Eq. 111|) 
describe two-particle inter-site transitions (^4,^3) «-» 
(k2, k\) of the renormalized fermions. For a given transi- 
tion, k is simply the number of virtual steps that have to 
be made by the original fermions. The steps are counted 
off from the configuration where two such fermions oc- 
cupy neighboring sites, and each step is a transition by 
one of the fermions to a nearest site. In other words, 
the original fermions go first from sites (k^k^) to sites 
(p,p+l) and then to {k%, k\) (we assume for concreteness 
that ks > &4 and k\ > k^)] the value of p is chosen so as 
to minimize the number of steps. 

To make the meaning of k even more intuitive we give 
examples of some >c = 4 transitions. For the initial and 
final states (n, n+l) and (n— 2, n+3) one of the sequences 
of steps of the original fermions is (n, n + 1) — > (n, n + 
2) — > (n— l,n+2) — * (n—l,n+3) — ► (n— 2, n+3), whereas 
for the initial and final states (n, n + 2) and (n — 1, n + 3) 
one of the sequences is (n,n + 2) — > (n,n + 1) — > (n — 
1, n + 1) — > (n — 1, n + 2) — > (n — 1, n + 3) [the energy 
denominators must be obtained directly from Eqs. JSJ, 

03]. 

It follows from the above argument that, for large it 
and a 3> ath, the transition matrix element 



Vfeife; 



fe 3 fc 4 



if* 



for 



x > 1. 



(14) 



The transitions of renormalized fermions are not limited 
to nearest neighbors. However, from Eq. 1)14(1 ■ the ampli- 
tudes of transitions over many sites are small and rapidly 
decrease with the number of involved virtual steps. The 
hopping integral for the renormalized fermions is VJA. 
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In higher orders of the perturbation theory, the inter- 
action leads also to many-particle transitions. The 
overall transition amplitude is determined by the total 
number of involved virtual single-particle steps. 

In order to localize many-particle excitations, one 
has to suppress combinational many-particle resonances 
keeping in mind that, for localization, the effective hop- 
ping integral must be smaller than the energy detuning 
of the initial and final on-site states. Because of the large 
number of possible resonances, we do not have an ana- 
lytical proof of many-particle localization for our energy 
sequence @. Instead we used numerical analysis, as de- 
scribed in the next section, which enabled us to demon- 
strate strong localization of stationary states in a chain 
of a limited size. 



B. The many-particle inverse participation ratio 

As mentioned earlier, a good indicator of strong local- 
ization, which applies to both single- and many-particle 
stationary states, is closeness of the IPR to one. In this 
section we present numerical data on the IPR obtained 
by diagonalizing the Hamiltonian (0) in the presence of 
several excitations. The Hamiltonian is a sparse matrix, 
which is separated into uncoupled blocks with different 
numbers of excitations. We focused on the block in the 
middle of the energy spectrum. For a given chain length 
L, it has L/2 excitations and therefore the largest number 
of states. This is the worst case, in terms of localization. 
Chains with L = 10, 12, and 14 were studied. The results 
were similar. We present the data for L = 12, in which 
case the total number of states is 924. 

The IPR as a function of the dimcnsionless parameter 
a of the on-site energy sequence (@J is shown in Fig. 0] 
The results refer to two values of the dimensionless ratio 
of the hopping integral to the interband distance J/h 
and several values of the dimensionless parameter of the 
particle interaction A. 

We start the analysis with the region a — > 0, where 
the on-site energies alternate between ±h/2. For 
large h/ J the single-particle energies form bands of width 
J 2 /2h. In the neglect of mixing of these bands, the many- 
particle wave functions can be found using the Bethe 
ansatz. For A = the many-particle energy spectrum 
consists of bands that are determined by the number of 
particles in each of the two single-particle bands. Be- 
cause the states are of the plane-wave type, the IPR is 
large, with (7 6 ) « 78 for h/J = 20 and L = 12, see Fig.H 

For a — > the IPR decreases with the increasing pa- 
rameter of the particle interaction A. This happens, 
because for large A ^> J/h (but A <C h/J) the en- 
ergy bands at ±h/2 split into subbands depending on 
the number of particle pairs, triples, etc. on neighboring 
sites: for example, the energy of a pair on neighboring 
sites differs from the energy of a dissociated pair by J A. 
Since the number of states in a subband is smaller than in 
the whole band, such splitting reduces the average IPR. 



100 




FIG. 4: (color) The IPR for 6 excitations on the first 12 
sites of the chain • The reduced bandwidth of the energy 
spectrum is h/J = 20 (top panel) and h/J = 10 (lower panel). 
The purple, red, green, and blue curves refer to the coupling 
parameter A = 0,0.3,1, and 3, respectively. The peaks of 
(Ie) for A = are a single-particle boundary effect. The 
insets show the maximal Iq for A/0. Sharp isolated peaks 
of ^6 max vs q result from the hybridization of resonating on- 
site many-particle states. 



The decrease of the average IPR with increasing A for 
a — > and A <c h/J is seen in Fig. 21 We note that the 
mere separation (by h) of the single-particle energies of 
neighboring states is not sufficient for strong localization. 
Inside each subband, the number of resonating states in 
a long chain is still very large. Localization requires that 
not only nearest neighbor single-particle energies, but 
also energies of remote sites be tuned away from each 
other. This should happen for sufficiently large values of 
the parameter a in Eq. 

This argument is confirmed by the data in Fig. ^ For 
given A, the IPR decreases as a whole with increasing 
a in the region where the single-particle bands are well 
separated, a < 0.4. This is the expected consequence of 
elimination of energy resonances. 

In the region 0.2 < a < 0.4, except for narrow peaks, 
we have l6max ~ 1-09 for h/J = 10 and l6max ~ 1-02 
for h/J — 20, for A < 1. The values of (Iq) are even 
smaller, 1.04 and 1.01, respectively. This indicates that, 
in this parameter range, all states are strongly localized. 
For A = 3 and h/J = 20 we also have /e m „ w 1.01 away 
from the peaks; however, for h/ J — 10 the IPR becomes 
larger due to the many-particle resonances, which are 
discussed below (see also Appendix C). 

A distinctive feature of the many-particle IPR as func- 
tion of a is the onset of multiple resonant peaks, which 
can be seen in Fig. 0] They indicate that at least some 
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of the stationary states are no longer confined to a single 
set of sites, with the number of sites equal to the num- 
ber of particles. The peaks are due to hybridization of 
resonating on-site states. The hybridization occurs when 
the matrix elements of inter-site transitions in Eqs. 
(fT2*)> exceed the energy difference of the states. 

For 0.2 < a < 0.4, i.e., in the region of strong local- 
ization, and for h/J = 20 and chosen A < 3 we found 
that only two on-site states could become strongly hy- 
bridized in the section of the chain with 1 < n < 12. 
Hybridization of a larger number of states was weaker. 
A consequence of hybridization of utmost two states is 
that i6max % 2. For h/J = 10 and A = 3 the interstate 
coupling iff^l is stronger, and as a result three states can 
be strongly hybridized and a few more can be weakly 
admixed, leading to -?6max ~ 3 at resonant a. 

Because the interaction is two-particle, the strongest 
peaks of / ma x come from resonances between on-site en- 
ergies of two particles. They occur when the energy dif- 
ference 

Se = \e kl + £/t 2 - £fc 3 - £fc 4 | (15) 

is close to MJA with M = 0, 1, 2. Strictly speaking, we 
should use exact single-particle energies e' n © instead of 
e n in Eq. i|15|) . but the difference between these energies 
is small, see Eq. (110(1 . and it leads to a small shift of the 
positions of the resonances as functions of a. An expla- 
nation of the positions and widths of the narrow peaks 
of (Is) and l6 max seen in Fig. 0] is given in Appendix C. 

C. Broad-band two-particle resonances 

A special role in the problem of many-body localization 
is played by two-particle resonances that are not selective 
in a, i.e., exist in a broad range of a. For these broad- 
band resonances, the total energy difference between the 
initial and final on-site states is small, Se ^ J even for 
small J/h. They emerge already in the second order in 
J/h, when only two single-particle steps to neighboring 
sites are required, k = 2. The resonating on-site states 
are pairs (n, n + 1) and (n — 1, n + 2), i.e., 

£n + £n+l ~ £ri-l + £n+2- (16) 

If n and n + 2 are prime numbers, the energy difference 
Se = \e n + e n+ i ~ e n -i - £ n +2\ ~ a n ^h is extremely 
small for large n. 

More generally, the resonance l(16|) occurs for all n = 
6k — 1 with integer k. In this case 5e/h oc with £ > 4. 
Such Se is "anomalously small" for x = 2. The hopping 
integral J 3 A/h 2 becomes larger than Se even when we 
are already deep in the single-particle localization region, 
a > a t h- 

As a result of renormalization e n — » e' n , the energy 
difference Se — \e n + e n +i — £«— l — £n+2| is changed. 
For n — 6k — 1 to leading-order in J/h the renormaliza- 
tion of Se does not exceed ~ a 3 (J 2 /h) or ~ a 4 (J 2 /h) 



for odd and even k, respectively (it may also be propor- 
tional to higher power of a). The renormalized Se can 
still be much smaller than J 3 A/h 2 , and then the broad- 
band resonance persists. More many-particle broad-band 
resonances emerge for larger x. 

For the section of the chain with sites 1 < n < 12 the 
IPR is not much affected by the broad-band resonances, 
because even where n and n + 2 are prime numbers (5 
and 7) , they are not very large and the energy difference 
Hlfty is not exceedingly small. 

The situation changes if we consider other sections of 
the chain (0J of the same length, no < n < tiq + 11 with 
different no- Here we found that the resonances increase 
(Iq) up to 1.15 between the peaks, for h/J = 20, A = 1, 
and for all a in the interval 0.2 < a < 0.4. 

For a finite chain, the resonances can be eliminated 
order by order in k by shifting the energies of the ap- 
propriate qubits. Simple systematic modifications of the 
energies that work for k < 5 are discussed in the next 
section, see Eqs. (tT§l) . Both modifications bring the 

IPR back to smaller values. For example, in all sections 
of the chain that we studied the made (Iq) and l6 max 
equal to w 1.01 and w 1.02, respectively, which are the 
values we had for the section 1 < n < 12. This indicates 
that the localization of stationary states for the modified 
energy sequences is very strong. 

V. LIFETIME OF STRONGLY LOCALIZED 
STATES 

The problem of strong localization can be viewed also 
from a different perspective. In the context of quantum 
computing, it suggests a more appropriate formulation 
then the one based on the analysis of stationary states. It 
is also relevant for condensed-matter systems at nonzero 
temperatures. 

First we note that excitations in quantum computers 
and in condensed-matter systems have a finite coherence 
time i C oh- For QC's, this time has to be compared with 
the duration of a single- or two-qubit operation and mea- 
surement. The duration of a two-qubit operation is of 
order of the time it takes to resonantly transfer an exci- 
tation between the qubits, which is ~ J -1 . A single-qubit 
operation is often faster; however, the measurement can 
sometimes be slower. In most proposed realizations of a 
QC the coherence time exceeds the gate operation time 
by a factor < 10 5 . 

We define the localization lifetime t\ oc as the time it 
takes for excitations to leave occupied sites. Localization 
of excitations is only relevant on times ~ £ C oh- Then 
to have strong localization it suffices that ii oc > i co h- It 
follows from the estimate for i co h that the latter condition 
is met if 

tioc > lO'j- 1 . (17) 

The condition l|17fl must be satisfied for all on-site many- 
excitation states. It is this condition that imposes a con- 
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straint on the form of the energy sequence e n in an infi- 
nite many-particle system. 

The time t\ oc is determined by hopping between res- 
onant on-site states. It occurs through virtual tran- 
sitions via nonresonant sites. For a two-particle reso- 
nant transition, the minimal number of the needed vir- 
tual steps is given by the parameter x (|13[1 . Then from 
Eqs. I|14|> the hopping integral for a resonant tran- 

sition (&4,£;3) <-> (fc2,&i) is 

JAV klk2k3ki ~ JAK" 

for x > 1. Here, K is defined by Eq. K ps J/2ah, 
and i-T -C 1 in the region a/ath 1. 

The localization lifetime t\ oc is determined by the re- 
ciprocal maximal hopping integral for resonating many- 
particle states. In the case of an energy sequence of the 
type up to a fairly high number of virtual steps (< 5), 
of interest are resonances between two-particle on-site 
states. This applies to systems with an arbitrary number 
of particles; only those transitions matter in which up 
to two particles change sites. Indeed, transitions where 
three particles change sites emerge in the second order 
in the two-particle Hamiltonian Q12JI. and simple 

counting shows that they involve already at least 5 vir- 
tual steps. 

For resonant two-particle transitions t\ oc ~ 
minfJAif^] -1 . Therefore it strongly depends on 
the minimal value of x for all pairs of resonating initial 
and final on-site states, x = x m i n . To have t\ oc J that 
exceeds a given value, we must have an appropriate x m i n . 
This means that we should eliminate resonances between 
all states connected by x < >c m i n virtual transitions. 

A two-particle transition with odd x = 1, 3, 5, ... in- 
volves a change of the total number of occupied sites 
with even n (and by the same token, with odd n, too). 
Therefore, for the sequence |@J with a <C 1, the energy 
change in such a transition is Se ~ h. In this paper we 
consider the case of weak to moderately strong coupling, 
when A < 1. Then Se significantly exceeds the change of 
the interaction energy JA, 2 J A for h <C J. As a result, 
resonant two-particle transitions may occur only for even 
x. 

We will modify the sequence (0} to eliminate reso- 
nances with x = 2 and x = 4. In these 
and 6, respectively, leading to the localization time t\ oc ~ 
J^K- 4 and > J^K- 6 . 



A. Eliminating second order many-particle 
resonances 

The potentially resonant transitions with x = 2 are 

(n, n + 1) <-+ (n, n + 1 ± 2), (n - 2, n + 1) 
(n,n+l) (n-l,n + 2). (18) 

In the transitions listed in the first line of this equation, 
one of the particles in the pair moves by two sites in one 



or the other direction, whereas for the transition shown 
on the second line both particles move by one site. 

The number of occupied nearest sites in the transitions 
(|18|l can change by one or remain unchanged. There- 
fore the maximal change of the interaction energy is JA. 
Second-order resonances will be eliminated if the detun- 
ing of the on-site energy differences Se for the transitions 
CB> is 

Se > JA. 

This means that we need a zero-energy gap of an appro- 
priate width in Se. 

We note that this is a sufficient, not the necessary con- 
dition. In principle, it would suffice to have narrow gaps 
at Se — and J A. These gaps should just be broader 
than the tunneling matrix element and than the energy 
shifts due to occupation of next nearest neighbors. For a 
specific finite-length section of a chain this may be more 
practical. However, here we are interested in an infinite 
chain, and we want to demonstrate that even for such a 
chain all resonances with x < 4 can be eliminated. 

To create the zero-energy gap, sequence Q has to be 
modified. The modification has to eliminate, in the first 
place, the "anomalous" broad-band resonances for transi- 
tions (n, n+1) <-> (n— 1, ra+2) with n = 6k — 1 discussed 
before. This will, of course, also eliminate resonances 
where n and n + 2 are prime numbers. A simple and 
sufficient modification is a constant shift of e n for each 
6th site, 

e" ld = e n + {h/2)a! for n = 6k, (19) 

while e" ld = e n for n ^ 6k. 

For the modified sequence (|19|) . the gap in the on- 
site energies for the 2nd-order transitions H18|) is Se ~ 
a 2 h,a'h/2 to leading order in a. A more accurate es- 
timate is minSe/h rj a 2 — a 3 ,a'/2. We assume that 
a 2 < a' <C a, in which case no new resonances are cre- 
ated for the transitions (|18|) as a result of the modification 

It follows from the above estimate that, for an infi- 
nite chain and an arbitrary number of particles, all res- 
onant transitions with x < 4 will be eliminated pro- 
vided JA/h < a 2 — a 3 ,a'/2, . Then the localization 
time ti oc ~ 10 5 J _1 already for h/J — 30, a = 0.3, 
a' w 0.1-0.2, and A < 1. 

The effect of the on-site energy modification (|19|l on 
localization of a many-particle state is seen from Fig. [S] 
This figure shows the dynamics of the system prepared 
initially in the state |$(416, 419, 420, 422, 423, 424)). In 
the case of the original sequence (@J , this state strongly 
hybridizes with the state |$(416, 418, 421, 422, 423, 424)) 
for all a of physical interest, a < 0.4. This happens 
because the difference of on-site energies £419 + £420 and 
£418 + £421 is Se ~ a 41& h. The hybridization results in 
oscillations of the amplitude of the state, as seen from 
Fig. For the modified sequence l|19f) the resonance 
is eliminated, and the amplitude remains constant over 
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FIG. 5: (color online) Time evolution of the squared ampli- 
tude \A\ 2 of the on-site state |$(416, 419, 420, 422, 423, 424)) 
in a 12-site section of the chain between n = 415 and n — 426. 
The oscillating line refers to the original sequence (3) with 
a = 0.25. The nearly constant line refers to the modified se- 
quence 11911 with a — 0.25, a' = 0.22. In both cases h/J = 20 
and A = 1. 



a time > 10 6 J~ 1 . This illustrates the onset of strong 
localization. We note that the localization time of this 
particular state turns out to be longer for the modified 
sequence (|19|l than the worst-case estimate given above. 



B. Eliminating 4th order resonances 

The localization time is further dramatically increased 
if k = 4 resonances are eliminated. The potentially res- 
onant 4th order transitions are 



(n, n + 1) «-> (n — 2, n + 3), 
(n,n + l) <-> (n + 2,n + 3), 
(n, n + 3) «-> (n — 1, n + 2), 
(n, n + 3) «-» (n — 2, n + 1), 



(20) 



and 



(n,n + 1) <-» (n- 1,71 + 4), 
(n, n + 1) <-» (n — 3, n + 2), 
(n, n + 1) <-> (n, n + 1 ± 4), (n ± 4,n + 1). (21) 

In the last line of Eq. i|21|) we list transitions where one 
of the particles in the pair moves by 4 sites, whereas in 
all other transitions both particles move away from their 
sites. 

For the modified sequence 119fl . the minimal energy 
change in the transitions (|20l) . I|21|) is minfc ~ a 3 h, to 
leading order in a, provided a' 3> a 3 . The value of a' has 
to be in such a range that the modification (|19fl does not 
lead to extra resonances between the on-site energies for 
the states 11201 and (I21|) . The "dangerous" combinations 
in Se/h are \a-a'/2\, \2a-a'/2\, \a 2 -a'/2\, \2a 2 -a'/2\, 



to leading order in a. We will choose a, a' so that all of 
them exceed min Se/h w a 3 . 

Fig.Elshows how the modification l|19|) leads to a zero- 
energy gap in Se. We plot fe™ d for all transitions lt2U|l . 
(|21|l with one of the particles on site n (with n > 2). 
Therefore we show all potentially resonant transitions 
with k < 5. 
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FIG. 6: (color online) The low-energy part of the two-particle 
energy differences Se n /h 1151 for all transitions with >c < 5 in 
which one of the involved particles is on the nth site (n > 2). 
The data refer to a = 0.25. The left panel corresponds to the 
sequence (|1J. The right panel refers to the modified sequence 
119t with a' = 0.22 and shows the zero-energy gap in Se. 



The left panel in Fig. H3 shows that, for the initial se- 
quence (@J), there is practically no gap in the values of 
5e at low energies. The right panel demonstrates that 
the correction l|19f) leads to a zero-frequency gap. The 
gap depends on the values of a and a' . For the specific 
parameter values in Fig. we have 6e/h > 0.01. 

We have checked numerically that the gap persists for 
a much longer chain than shown in Fig. [U with n from 
2 to 10,000. This is more than enough to prove that 
the results apply to an infinite chain. Indeed, in e™ d 
the terms oc ofl with different q are repeated periodically 
with period 2(q+ 1). Therefore a sequence of terms with 
a q for q < 6 is repeated periodically with period equal to 
twice the least common multiple of all q+ 1 < 7, which is 
2 x (3 x 4 x 5 x 7) = 840 (cf. Appendix A). For a = 0.25 
we have a 6 ~ 2 x 10~ 4 . The contribution to e" ld //i of 
all terms of higher degree in a is then < 2 x 10~ 4 . This 
means that, to accuracy better than 2 x 10~ 4 , the results 
on the gap for an infinite chain will coincide with the 
results for a chain of 840 sites. 



It follows from the discussion above that, for 2JA < 
O.Ol/i, all particles will remain localized on their sites 
for the time t\ oc ~ ( JA)~ 1 K ~ 6 /a [we have taken into 
account here that the hopping integral for transitions 
with x = 6 is limited by ~ JA( J/2h) 6 a~ 5 rather than 
JA( J/2h) 6 a~ 6 , as would be expected from the asymp- 
totic expression (|14|) ]. For h/J = 50 and a = 0.25 this 
gives an extremely long localization time, t\ oc J > 10 10 . 
However, this estimate requires that the coupling be 
weak, A < 0.25 for the used parameter values. 
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1, Extension to stronger coupling 

The previous result can be easily extended to larger A 
without increasing h. To do this we note that the change 
of the interaction energy in transitions i|20|) is actually 
limited to J A rather than 2 J A, which is the case for the 
transitions (|21[1 . Therefore the gap for the transitions 
(|21|l should be twice as large as for the transitions (|20l) . 
Two first transitions (|21l) have a gap > a 2 h to leading 
order in a, but for the last one minfc = o?h. This latter 
gap may be increased by choosing a somewhat different 
modification of the on-site energy spectrum. In contrast 
to Eq. I|19(l . we will shift the energy of each 3rd site, 



~md 



(ty4)/3[l + 3(-l)» 



for 



3fc, (22) 



while £™ d = e n for n ^ 3k. 

Eq. (22)) corresponds to shifting e^k up by f3h/2 or 
down by (3h depending on whether k is odd or even, re- 
spectively. The parameter (3 should be much larger than 
a 3 , to open a gap ~ min(a 2 /i, fih/2) in <5e™ d for the tran- 
sitions (|21|) . At the same time, it should be chosen so as 
to avoid creating new resonances, similar to the modifi- 
cation ljT3|) . 

It is straightforward to show that the modification 11' I'D 
leads to a zero-energy gap in Se. For a = 0.25, (3 = 0.1725 
the gap exceeds O.Ol/i for the transitions (|18fl . I|20|) with 
energy change < J A, whereas for the transitions H21|) 
with energy change up to 2JA it exceeds 0.02/i. This 
indicates that the results on the localization time t\ oc J > 
10 10 for h/J = 50 now apply for A < 0.5. 

As in the previous section, the extremely large localiza- 
tion time characterizes an infinite chain and an arbitrary 
number of interacting particles. We note that both mod- 
ifications of the original energy sequence, Eqs. (|19fl and 
(|22|l . are obtained analytically, by finding the leading- 
order terms in the energy differences for the transitions 
(|18fl . (|20() . and 1|21[) . The specific values of the parame- 
ters a, a', and /? are used just to illustrate the order of 
magnitude of the localization time. 



VI. STABILITY WITH RESPECT TO ERRORS 
IN ON-SITE ENERGIES 



The effect of errors on the gap can be modelled by 
adding a random term to on-site energies, i.e., replacing 



e™ d with 



-Dhr ri 



(23) 



Here, r n are random numbers uniformly distributed in 
the interval (—1,1), and D characterizes the error am- 
plitude. It should be compared with a s with different 
exponents s > 1. When D ~ a s it means that the ener- 
gies e n are well controlled up to terms ~ to leading 
order in a. 

From the above arguments it follows that, for a » 
a' > a 2 the gap should remain unchanged if D <C a 4 . 
This is because, for the modified energies e™ d , the terms 
~ a 4 drop out from the energy differences that we dis- 
cuss. For D ~ a 4 the gap should be somewhat reduced. 
For D ~ a 3 it should become significantly smaller than 
for D = 0, and it should ultimately disappear with in- 
creasing D. 

Numerical results on the gap Se as a function of log D 
are shown in Fig.Q The gap is calculated for two-particle 
transitions with >t < 5, as in Fig. |SJ In the lower panel 
the gap is scaled by its value in the absence of errors, 



i? = minfe° r 7mmfe™ 



md 



(24) 



The data refer to the same a, a' as in Fig. They are 
in full agreement with the above estimate. 
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In a real system, it will be impossible to implement 
sequence of on-site energies (QJ precisely. This is because 
these energies contain high powers of the small parameter 
a, while the precision to which they can be set and/or 
measured is limited. Therefore it is necessary to study 
localization in the presence of errors in e n and to find how 
large these errors can be before they cause derealization. 

We will address this problem by looking at the gap in 
the energy differences Se in the presence of errors in e™ d 
(|19fl . As long as this gap remains larger than 2JA for all 
resonant transitions with x < 5, the localization lifetime 
ti oc will remain large. The analysis can be immediately 
extended to the modified sequence (122[) as well. 



FIG. 7: (color online) Upper panels: all energy differences 
SeT/h = \eT + e c ™ - En" - eT 3 \/h for the transitions JTHJ, 
1201 . 12111 that correspond to the number of intermediate steps 
h < 5. The data refer to a = 0.25, a' = 0.22, and to a 
specific realization of the random numbers r n in Eq. \1'M . 
The boxes from left to right correspond to the values of the 
noise intensity D — a 3 in Eq. (IL'.'-il with s — 5,4, and 3. 
Lower panel: the scaled minimal gap R 12 11 as a function of 
the logarithm s — In D/ In a averaged over 10 realizations of 
noise. Error bars show the standard deviation of R. 

The results of Fig. [7| demonstrate that the localization 
persists even for relatively large errors in the on-site en- 
ergies. At least for the chosen a and a', errors in e n up 
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to ~ 0.4% (when D — a 4 ) lead to a change in the width 
of the energy gap by ~ 50%. 

The observed dependence on the noise strength sug- 
gests that, in the presence of noise, sequence (QJ, <(T§|) 
can be cut so that the terms oc a s with s > s C utoff are 
disregarded. The value of s cu toff depends on the noise, 
^cutoff = ha D I In a. As a result of the cutoff, the energies 
e n become polynomials in a of power < s cuto ff. From 
Eq. Q, these polynomials are periodic in n, with the 
period determined by twice the least common multiple 
of (2,3,..., s cu toff + !)■ For example, for s cuto ff = 5 the 
period in n is 120 (see also Appendix A). For such a long 
period and short-range hopping, excitations will stay on 
their sites for a long time compared to J . 



VII. ALTERNATIVE ENERGY SEQUENCES 

Neither the original sequence of on-site energies (QJ 
nor its modified versions Ijl9|) . I|22|) were optimized to 
maximize the IPR or the localization lifetime. Therefore 
it is important to compare them with other sufficiently 
simple sequences. This will be done in this section for 
two natural choices of e n . 



A. Period doubling cascade 

A simple way to move sites with close energies far away 
from each other is to make the energies form a "period 
doubling cascade" (PDC). It can be described by a one- 
parameter energy sequence; in what follows h is the en- 
ergy scale and a is the parameter, as in Eq. (QJ. 

In the PDC, the on-site energies are first split into two 
subbands, with nearest neighbors being in different sub- 
bands, but next nearest neighbors being in the same sub- 
band. The subbands differ in energy by a°h, to leading 
order in a [this is also the case in Eq. (0}] . Each of the 
subbands is then further split into two subbands of 4th 
neighbors. The leading term in the energy difference of 
these subbands is a^h. Each subband is then split again 
into two subbands of 8th neighbors. The leading term in 
their energy difference is a h. This period-doubling pro- 
cess is then continued indefinitely, for an infinite chain. 

The expression that describes the on-site energy se- 
quence for the PDC, e^ DC , can be conveniently written 
in terms of the coefficients jk(n) = 0, 1 of the expansions 
of site numbers n in base two, 

M(n)-1 
»= X] jk(n)2 k . 

k=Q 

Here, M(n) = 1 + Ll°g2 n \ is the number of integer digits 
of n in base 2. 
We set 



The energies e^ DC are shown in Fig. [5] It is instructive to 
compare this figure with Fig. ^ for the energies J3J. The 
overall band structure is similar, but the energy distribu- 
tion is much more regular for the PDC. For example, it 
is seen from the panels for 2000 sites that the minibands 
for the PDC sequence have approximately equal num- 
bers of states, in contrast to Fig.^ However, as we show 
later, the symmetry of the PDC is actually bad from the 
viewpoint of many-particle localization. 
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FIG. 8: (color online) The on-site energies £^ DC //i for 
a = 0.3. The left panel shows the energies for the sites 
n — 1,2,..., 50. States with close on-site energies are spa- 
tially separated. The right panel shows e^ DC //i for a much 
longer array, n = 1, . . . , 2000. The energy spectrum displays 
a multisubband structure. Because of the symmetry of se- 
quence 1251 1. the number of points in each subband is approx- 
imately the same. 

Spatial separation of sites with close energies in the 
PDC leads to effective localization of single-particle sta- 
tionary states. As in the case of sequence (@J , the values 
of (ii) — 1 are ~ 3 x 10~ 3 for h/J~20 near the minimum 
of (7i) over a, in agreement with the estimate (J7J). The 
minimum of (ii) is located at a « 0.1. 

The situation is different for many-particle localiza- 
tion. Here the high symmetry of the PDC leads to multi- 
ple many-particle resonances. For example, two-particle 
states (n, n + 1) and (n — 1, n + 2), which are coupled 
in second order in J/h, have equal energies whenever n 
is odd. The states [n,n + 1) and (n — 2,n + 3), that 
are coupled in the fourth order in J/h, have the same 
energies when n = 7 + Ak with integer k. We note that 
there were no exactly degenerate states for sequence @ , 
and the portion of two-particle states with close energies 
was much smaller. Therefore it is more complicated to 
find a correction to sequence (|25|l that would eliminate 
many-particle resonances. As a result, unexpectedly, this 
symmetric sequence is less convenient from the point of 
view of strong localization. 



B. Random on-site energies 

The case opposite to the highly symmetric sequence 
(|25|l is when the on-site energies e n are completely ran- 
dom. It is well-known that such randomness leads to 
single-particle localization of stationary states in a ID 
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chain. However, it does not lead to strong on-site local- 
ization of all states, because there is always a nonzero 
probability to have neighboring sites with energies that 
differ by less than J and therefore are hybridized. As ex- 
plained in the Introduction, hybridization is even more 
likely to happen in the case of many-particle states, be- 
cause it is more likely to have neighboring nearly resonant 
sites. 

A simple random sequence of on-site energies has the 
form 

el = Wr' n , (26) 

where r' n with different n are independent random num- 
bers uniformly distributed in the interval (0, 1), and W is 
the bandwidth. The results on the inverse participation 
ratio for sequence (|26|l are shown in Fig. [5] 
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FIG. 9: (color online) The normalized distributions of the 
IPR P(I) for the random energy sequence (12611 . The upper 
and lower panels refer to the many- and single-particle IPR's: 
6 excitations with A = 1 in the chain of length L = 12 and 
one excitation in the chain of L = 300, with 2000 and 80 
realizations, respectively. The left and right panels refer to the 
overall bandwidth of the on-site energy spectrum W/J = 13 
and 26. The narrow columns at I — 1 showed with dotted 
lines refer to sequence @, with h/ J = 10 and 20 for the left 
and right panels, respectively. The values of a in Eq. I^J were 
chosen so as to have the same bandwidths of on-site energies 
as the corresponding random sequences: in the lower panels 
a ~ (W — h)/W « 0.231 (as in an infinite chain), whereas 
in the upper panels a — 0.274 to allow for a comparatively 
small number of sites. The insets show log P vs. I. 

Localization of stationary states for the random se- 
quence (|26[1 can be characterized by the probability dis- 
tribution of the IPR P(I). This distribution was ob- 
tained by numerically diagonalizing the Hamiltonian @ 
for different realizations of the on-site energies l|26[) . 

As seen from Fig. both single- and many-particle 
IPR distributions display peaks near 1=1 for random 
on-site energies. This indicates that, for the wide energy 



bandwidths used in the calculations, most of the station- 
ary states of the system are strongly localized. How- 
ever, the distributions are broad and slowly decay on the 
tails. This means that many on-site states are strongly 
hybridized, that is the stationary wave functions spread 
over several sites. This is a consequence of multiple res- 
onances. The insets in Fig. [5] show that, at least for not 
too large /, the tails of P{I) are non-exponential. 

The typical distribution width decreases with the in- 
creasing bandwidth W of on-site energies. As expected, 
the distributions of many-particle IPR's (the upper pan- 
els in Fig. are much broader and their peaks near 1=1 
are much smaller than for the single-particle IPR's (the 
lower panels). 

The IPR distributions for the random sequence l|26|) 
differ dramatically from the distributions for the regular 
sequence (@J). The latter are narrow and concentrate in 
a small region of / close to I = 1, for the chosen band- 
widths, both for the single- and many-particle states, cf. 
Figs. |21 El This is another indication of strong single- 
and many-particle localization for sequence 



VIII. CONCLUSIONS 

In this paper we have explored two aspects of the prob- 
lem of strong many-particle localization. One is localiza- 
tion of stationary states. For one-particle states, it has 
been studied analytically. We found that the wave func- 
tions decay quasi-exponentially and obtained the bounds 
on the decay length. The numerical results on strong 
localization are in agreement with the theory. 

For many-particle stationary states, the localization 
has been analyzed numerically. Such analysis is unavoid- 
ably limited to small chains. For a 12-site chain we have 
found that, for the on-site energy sequence l|19|) with the 
ratio of the single-site energy bandwidth to the hopping 
integral » 28 [h/J = 20,a = 0.25, a' = 0.22 in Eq. <pi) ]. 
the inverse participation ratio differs from its value for 
the case of fully localized states by < 2%. This devia- 
tion is due to a small nonresonant admixture of the wave 
functions of neighboring sites. 

A different approach is based on studying the lifetime 
of on-site states. It is sufficient to have a localization 
lifetime ti oc that exceeds the coherence time of the exci- 
tations. We have shown that such t\ oc can be achieved in 
a chain of an arbitrary length and with an arbitrary num- 
ber of excitations. For the explicit construction of on-site 
energies <|19[) . resonant transitions that lead to delocal- 
ization require at least 4 virtual nonresonant steps, and 
for comparatively weak particle-particle coupling even 6 
virtual transitions. Even in the 4-step case this leads to 
the ratio of the delocalization rate to the inter-site hop- 
ping integral ~ 10~ 5 already for the decay length of the 
amplitude of virtual transitions ~ 0.35, or for K ~ 0.06. 

The results on localization can be tested with artifi- 
cial condensed-matter structures, as it was done in the 
studies of the effects of quasiperiodicity [2(j (see also 
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Refs. ll8L ll QT> . A discussion of experimental implementa- 
tions is beyond the scope of this paper. Good examples 
are proposed models of a quantum computer; we note, 
however, that studying localization does not require op- 
erations on qubits, and therefore does not require a fully 
operational quantum computer. 

In terms of quantum computing, an advantageous fea- 
ture of sequence (J2J and its modification (|19fl is that one 
radiation frequency can be used to resonantly excite dif- 
ferent qubits. This can be achieved by selectively tun- 
ing them to this frequency without bringing neighboring 
qubits in resonance with each other, or by sweeping the 
frequency of the targeted qubit through the radiation fre- 
quency and having a Landau-Zener-type interstate tran- 
sition. A two-qubit gate can be conveniently done by 
selectively tuning neighboring qubits in resonance with 
each other and having a Landau-Zener type excitation 
swap 16]. 

In many cases localization is also a prerequisite for a 
projective measurement. This happens when the mea- 
sured quantity is the probability for each qubit to be in 
the excited state. Often a measurement is much slower 
than the time J -1 of resonant hopping to a nearest site; 
then <i oc should exceed the measurement time. In our ap- 
proach, localization does not require refocusing 4|, which 
is not always easy to implement and which is sometimes 
incompatible with slow measurement. 

Our energy sequence does not lead to derealization 
of stationary states due to indirect resonant n «-> n + 2 
transitions. Such transitions undermine the approach to 
quantum computing with "always-on qubit interaction" 
proposed in Ref. Il2 . where next nearest neighbors are 
tuned in resonance. In terms of the localization lifetime, 
the width of the energy distribution 1)19(1 is orders of mag- 
nitude smaller than the width that would give the same 
tioc in the approach 12]. There the localization lifetime 
would be ~ h/J 2 , where h is the energy difference be- 
tween nearest neighbors. This difference must exceed J 
by a factor ~ 10 5 in order to have the same localization 
lifetime as in our approach. Compared to Ref. |Tlj , a po- 
tential advantage of our approach is that the interaction 
does not have to be ever turned off, and no multi-qubit 
encoding is necessary for operating a QC. 

The presented scheme can be extended to systems with 
long-range coup ling . This is another distinction from 
the approaches [11L Il2| that substantially rely on near- 
est neighbor coupling. For several proposed QC's the 
interqubit coupling is dipolar for a few near neighbors 
and becomes quadrupolar or falls down even faster for 
remote neighbors [a, |6j, |9( . Long-range interaction makes 
transitions over several sites more probable. We leave de- 
tailed analysis for a separate paper. Here we note that, 
for sequence l|22|) and for single-particle transitions, not 
only hopping over 2 or 4 sites is not resonant and does 
not lead to derealization, but even hopping over 6, 8, 
or 10 sites is nonresonant as well. For all these transi- 
tions, the energy difference is at least ~ a 2 h or ~ f3h/4 
[for sites separated by an odd number of positions, the 



energy difference is always large, ~ h]. 

Our results provide proof of principle of strong on-site 
localization of all states, independent of the system size 
and the number of particles. We have not addressed the 
question of optimization of the energy sequence, so that 
maximal localization lifetime could be obtained for a min- 
imal bandwidth of on-site energies. For a finite-length 
chain the optimization problem can be approached using 
Eq. J3J as an initial approximation and adjusting energies 
of several specific sites. 

In conclusion, we have proposed a sequence of on-site 
energies (@J and its modifications (|19|) , l|22|) that result in 
on-site localization of all single- and many-particle states 
of interacting spins or fermions in an infinite chain. The 
sequence Q has low symmetry, which allows eliminat- 
ing resonances between the states to a high order in the 
hopping integral. In turn, this leads to a long lifetime of 
localized many-particle states. When second-order reso- 
nances are eliminated, the lifetime exceeds the reciprocal 
hopping integral by 5 orders of magnitude provided the 
bandwidth of on-site energies is larger than the inter- 
site hopping integral by a factor of ~ 40. We show that 
it can be further significantly increased by eliminating 
fourth-order resonances. The proposed energy sequence 
is stable with respect to errors. The results apply to 
scalable quantum computers with perpetually coupled 
qubits. They show that, by tuning qubit energies, exci- 
tations can be prevented from delocalizing between gate 
operations. 
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APPENDIX A: EXPONENTIAL DECAY OF THE 
TRANSITION AMPLITUDE 

In this Appendix we give a rigorous proof of the quasi- 
exponential decay of the amplitude K n (m) (JSJ of the 
transition from site n to site n + m and establish bounds 
on the decay length. We show that, for sequence (0J, in 
the limit of small a and for m — > oo 

(X~ v l m < K n (m){2h/J) m < a~ v u m . (Al) 

We find that v L > 0.89 and v u < 1.19. 

In order to simplify notations we introduce dimension- 
less energies e n (a) — 2s n /h. From Eq. (J3J 

n+l 

e n (a) = (-1)" - £(-l)L»/kJ a *-i. (A2) 

fe=2 
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We also set J/2h = 1. Then K n {m) = l/\Q n (m)\, where 



Qn (m) = ]J I £ "+- (a)-e n (a)}. (A3) 

s=l 

From Eq. l)A2j) . Q n (m) is a polynomial in a. 

For a polynomial -P(a) we define by lowdegP(a) the 
multiplicity of the root a — 0, i.e., the lowest power of 
a in the polynomial. The exponent v that characterizes 
the decay of K n (m) © is given by v — lowdeg Q n (m)/m 
for m — > oo. 

The data of the numerical experiments presented in 
Fig. 0show that 0.894 < lowdeg Q n (m)/m < 1.12 inde- 
pendent of n for large m. 

To obtain an analytical estimate we rewrite Eq. (|A3|) 

as 

m 

lowdeg Q n (to) = lowdeg [e n+s (a) - e n (a)]. (A4) 

s=l 

Each term lowdeg [e„ +s (a) — e n (a)] is an integer between 
and n + 1 . 

In order to find bounds for lowdeg Q n {m) we will es- 
timate how many terms lowdeg [s n + s (a) — £n(a)] exceed 
a given i, < i < n + 1. For each i we have a subset 
S n m{i) of the values s that satisfy this condition, 

S nm (i) = {s\l < s < to, lowdeg [e n+s (a) - e„(a)] > i}. 

(A5) 

The number of elements in S nm (i) is denoted by h nm (i). 
This is the number of polynomials e n + s (a) — £„(a) whose 
expansion in a starts with a k with k > i. In what follows 
for brevity we drop the subscripts n, m and use S(i) and 
h(i) for S nm (i) and h nrn (i). 

It follows from the definition that 

Mo) > M 1 ) > H 2 ) > ■ • • > 

By construction 

n 

lowdegQ n (m) = 2_. i[h(i ~ 1) — h(i)} 

+(n + l)fc(n) = X£ =0 fc(i) (A6) 

From Eq. (|A6|) we see that the upper and lower bounds 
on lowdeg Q n (to) are given by the sums of the upper and 
lower bounds of h(i). 

In what follows we will use the standard notations: 
liminf (limsup) means the lower (upper) limit of a se- 
quence, LCM{ix 5 • • • j v} is t ne least common multiple 
of integers i%, . . . ,i r , and GCD{ii, . . . , i r } is the great- 
est common divisor of ii, . . . , i r . We will also denote by 
[e(a)]fc the coefficient at a k in the polynomial e(a), i.e., 

n 

£n{a) = ^[e„(a)] fc a fc . (A7) 

fc=0 



1. Lower bound 

In this section we obtain the lower bound of 
lowdeg Q n (to). The main statement is the following 
lemma. 

Lemma A.l The lower bound has the form 
liminf 1OwdegQ " (m) >0.89. 

m^oo TO 

Proof. Consider first the constant term [e n (a)]o in 
Eq. l|A7p . By definition, [e„(o;)]o = [e„ + 2(a)]o, and 

1 A I ( \ ( \\ J fOT ° dd S ' 

lowdeg((T„ +s (a;) - e n {a)) = < 

> 1 tor even s. 

Hence, we immediately obtain a simple lower bound 
lowdeg Q n (to) > \m/2\ for large to (in what follows we 
always imply to — > oo). 

We compute h(0), h(l), etc, using that the coefficients 
[e n (a)]i are periodic in n with period 2(i + 1). Indeed, 
from Eq. JX§J|, 

[e n (a)], - (_i)L«/«+i)J 

= (-l)L[n+2(i+i)]/(i+i)J =[eTH . 2(<+1)(a)]i . 

Therefore, the sets of coefficients 

{[e„(a)] , [e„(a)]i, . . . , [e„(a)]i} are also periodic in 
n, but with the period T 4 = 2LCM{2,3, ...,« + 1}. This 
is illustrated by the table 

£\{a) = —1 —a 

e 2 (o;) = 1 +a —a 2 

£3 (a) = -1 +a +a 2 -a 3 

64(a) = 1 —a +a 2 +a 3 —a 4 

85(a) = —1 —a +a 2 +a 3 +a 4 —a 5 

e 6 (a) = 1 +a -a 2 +a 3 +a 4 +a 5 ~a e 

e 7 (a) = -1 +a -a 2 +a 3 +a 4 +a 5 +a 6 -a 7 

In order to estimate h(i) we need two technical state- 
ments. 

Lemma A. 2 Let ao,k,T be any integers such that 2k 
does not divide T. Consider any 2k/ GCD{T, 2k} 
consecutive elements of an arithmetic progression cij — 
oo + j ■ T, and set bj — [oj/k\ mod 2. 

Then, at least [fc/ GCD{T, 2fc}J integers bj are egual 
to 0, and at least the same number of bj are equal to 1. 

Proof. Since 2k does not divide T, sequence aj mod 2k 
is cyclic in the interval [0,2fc — 1]. This sequence con- 
tains exactly 2k/ GCD{T, 2k} distinct elements. On av- 
erage, half of them (at least \p/ GCD{T, 2k}\) are less 
than p, and another half (the same number) are larger 
or equal than p. This means that there are at least 
[k/ GCD{T, 2k}\ integers bj that are equal to and at 
least the same number of bj that are equal to 1. Q.E.D. 

The next statement is a corollary of lemma IA.2I and 
we skip the proof. 
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Corollary A. 3 Let 2k does not divide T. Consider p 
consecutive elements of the arithmetic progression aj — 
a o + j ' T and set bj — \aj /k\ mod 2. 

Then at least [k/ GCD{T, 2k}\ [p ■ GCD{T, 2fc}/2fcJ 
integers bj are equal to 0, and at least the same num- 
ber of bj are equal to 1 . 

We are now in a position to finish the proof of 
Lemma IA.1I We notice first that, for n — dj in the 
expression l)A2[l . the coefficient bj for given k determines 
the sign of the term a fc_1 in e n (a), that is (— l) bj — 
[e n (a)]k-i- The number h(i) gives the probability that, 
for all k < i, the polynomial £ n + s (a) has the same bj as 
e n (a). 

We will now estimate h(i) with i — 1, . . . , 4 and start 
with h(l). We note that s e 5(1) if and only if s E 5(0) 
and [e„_|_ s ]i = [e n ]i< The second condition means that 
[s/2\ mod 2 = 0. By construction (|A5|) . for m — > oo 
the set 5(0) is formed by all numbers s of the same 
parity as n. This means that 5(0) is an arithmetic 
progression with period To = 2. We take p consecu- 
tive elements si, . . .s p of it and use Corollary IA.3I with 
k = 2, because we are interested in the coefficients 
[e n+Si ]i = (-l)L(»+^)/2J in By Corollary FOl 

since To is not divisible by 2k, for at least [p/2\ subscripts 
Si the coefficients [e n+Si (a)]i = 1, and [e„ +Si (a)]i = — 1 
for at least the same amount of subscripts Sj, i.e., approx- 
imately half of the coefficients [e n + s< ( a )]i coincide with 
[e n (a)]i. Hence, h(l) > h(Q)/2 as to — > oo. Substituting 
h(0) = m/2 we obtain h(l) = m/4. 

Similar arguments can be applied to estimate h(2). 
This requires finding a portion of the set 5(1) which 
forms S(2). The set 5(1) is a nonempty disjoint union 
of arithmetic progressions with period T\ — 4. We 
will apply Corollary IA.3I with k — 3 to each of these 
progressions and use p = h(l). This gives h(2) > 
L3/2J L(l/3)(m/4)J = [m/12\, or h(2)/m > 1/12 for 

TO — > OO. 

In the same way we obtain h(3)/m > 1/24 and 
h(4)/m > 1/60 as m — » oo. Therefore 

lowdeg Q n (to)/to >i + I + l + l + i-> 0.89, 

which finishes the proof of the lower bound. 



2. Upper bound 

We start with the proof of the following rough esti- 
mate: 

Lemma A. 4 An upper bound has the form 

lowdegQ„(m) 22 
hmsup < — < 1.47. 

n — *oo rn 15 



Proof. Taking in the rhs of Eq. <| A6|> the sum from to 
oo we obtain 



To find an upper bound on h(i) we will use the follow- 
ing consequence of Lemma IA.2I 

Corollary A. 5 Let 2k does not divide T . Consider p 
consecutive elements of the arithmetic progression aj — 
oq + j ■ T and set bj = \aj/k\ mod 2. 
Then at most 



2k 



(oc I • / J 



GCD{T, 2k} 



p-GCD{T,2fc} 
2k 



1 



lowdeg Q n (m) < h(i). 



(A8) 



i=0 



integers bj are equal to 0, and at most the same number 
of bj are equal to 1 . 

Using the same arguments as before, by corollary I A. 51 
we obtain for m — > oo the following upper bounds for 
h(i) 



h(0) < m/2, 

h(l) < m/4, 

h(2) < m/6, 

h(3) < m/12, 

h{4) < m/20. (A9) 

Recall that h(4) > h(5) > h(6) > h(7). Similarly, for 
q > 2 we have 

h{2 q )>h(2" + l)>...h(2 q+1 -l). 

Therefore we can replace the terms 

h(2i + l),...,/i(2« +1 - 1) in the rhs of Eq. (HJ by 
h(2 q ), which leads to the following upper bound for 
lowdeg Q n (m), 

oo 

lowdeg Q n (m) < h(Q) + h(l) + h(2) + h(3) + ^ 2 q h(2 c '). 

q=2 

(A10) 

This reduces the calculation to finding upper bounds on 
h(2 q ). 

We will now obtain a recurrence relation for h(2 q ). 
First we notice that, for K + 1 being a prime number, we 
have from Corollary IA. 51 

h{K + l) - 2{K+l) h[K) - (AU) 

For all primes K + l > 7 we have (K +2)/2(K +1) < 4/7. 
Therefore h{K + 1) < jh(K) for K > 6. We also note 
that, for all positive integral q, 

h(2 q )<±h(2 q -l). 

Now we recall the distribution law for primes in the 
intervals. The following statement is called Bertrand's 
postulate (or Tchebychev' theorem) (see Ref. 152): 
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Theorem A. 6 There is at least one prime between M 
and 2M for any positive integer M. If M > 3, there is 
always at least one prime between M and 2M — 2. 

In particular, there is at least one prime between 2 q 
and 2 q+1 — 1 for any positive integer q > 2. With this 
statement, taking into account the previous estimates, 
we obtain 




for all positive integral q, or 

M2" +2 )< 0)"^. (A12) 



Substituting inequalities HA9jl. <A12(1 into (|A10|) . we 
obtain 




for to — > oo. 



This gives 

22 

lowdeg Q n (to) < — to < 1.47m (A13) 
15 

The last inequality is an explicit asymptotic upper 
bound. Q.E.D. 

We now provide a sharper upper bound. We will use 
the same method, but instead of Tchebyshev' theorem 
we will apply Erdos theorem. 

Lemma A. 7 (Sharper bound) 

lowdeg Q n (to) 
limsup < 1.19 



Theorem A. 8 (Erdos) There exist at least one prime of 
the form 4fc + 1 and at least one prime of the form 4fc + 3 
between M and 2M for all to > 6. 

For all primes that exceed 16 we have in Eq. (|A11|) 
{K + 2)/2{K+ 1) < 7/13. Therefore, by repr oducing the 
arguments that led to the inequality l|A12l) . but using 
now the relation 1A11|) twice based on the theorem IA. 81 
we obtain 



h(2 q+i ) < 



h(16). 



(A15) 



Proof. Following the same pattern as in Lemma I A. 41 
above we extend the explicit list of inequalities (|A9|) . 
As to — > oo, 



h(5)/m 


< 


1/20, 


h(6)/m 


< 


1/35, 


h(7)/m 


< 


1/70, 


h{8)/m 


< 


1/126, 


h(9)/m 


< 


1/126, 


h(10)/m 


< 


1/231, 


ft(ll)/m 


< 


1/231, 


h(12)/m 


< 


1/429, 


h(13)/m 


< 


1/429, 


h(U)/m 


< 


1/429, 


h(l5)/m 


< 


1/858, 


h(16)/m 


< 


3/4862 



(A14) 

To obtain a sharper upper bound we recall the follow- 
ing result by Erdos [23J. 



Sub stituting inequalities (|A9I) . (|A14|I . (|A15I into 
(|A10(1 (where now the terms up to h(15) are taken into 
account explicitly, and the sum runs from q = 4) we ob- 
tain 



lim sup 



lowdeg Q n . 




Evaluating the rhs, we obtain 



lowdeg Qn (to) 
limsup < 1.19. 



(A16) 



Q.E.D. 
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APPENDIX B: THE INTERRELATION 
BETWEEN THE ENERGY SPECTRUM 
PARAMETERS a AND J/h FOR FIXED IPR 

In this Appendix we outline another way of looking 
at the effect of the band structure of sequence 011 
localization. It applies to single-particle stationary states 
and is based on varying h/J and finding such energy 
spectrum parameter a that would keep the IPR constant, 
i.e., 

(h) = (h(a, J/h)) = const. (Bl) 

The average IPR is large, w L/3, when the spread of 
the on-site energies ah is small compared to the hopping- 
induced bandwidth J 2 /2h of the bands at ±h/2. When 
ah becomes comparable to J 2 /h, a part of the states 
become localized with localization length smaller than 
the chain size, but still there remain states of size ~ L. 
For such states I\\ oc L. Their portion depends on 
ah/(J 2 /h). Therefore one may expect that, for large 
and for a given chain length, a should vary with 
J/h as (J/h) 2 . 

Another scaling region of a(J/h) as given by Eq. fBTJl 
may be expected to emerge for a close to the thresh- 
old value, ath < a <C 1, but far away from the strong- 
localization range of a, where (Ii) — 1 ~ J 2 /h 2 [cf. 
Eq.0]. 

For a close to ath, the wave functions have compara- 
tively small-amplitude tails that spread over a long dis- 
tance and are nearly exponential at large distances, as 
given by Eq. J^J. If the decay were purely exponen- 
tial, i.e., the tail of the wave function centered on site 
n were of the form ipn+m = K^ m 'tp n , we would have 
(Ji) - 1 re 4\K\ 2 for K Pa J/2ah < 1. From jBlJ, this 
condition gives scaling a oc J/h. The nonexponential de- 
cay of the wave functions at small to moderate distances 
(numerically, for |m| ~ 4 — 8) leads to deviations from 
this simple scaling. 

Numerical results on the dependence of a on J/h as 
given by Eq. (|B1(1 are shown in Fig. ^3 The data for 
(ii) ~ L/3 show the expected scaling a oc (J/h) 2 . On 
the other hand, in the range (Ji) - 1 « 0.1 — 1 the value 
of a scales as J/h. This scaling applies only for a > ath, 
i.e., for ah/ J > 1/2. The value of ah/ J as given by 
Eq. I|B1I) increases with decreasing (ii). 

We note that, for large h/J^ 100 and small (Ii) — 1, 
the IPR (Ji) as a function of a displays small oscillations. 
This leads to multivaluedness of the roots a of the equa- 
tion (Ji) = const. The roots are numerically very close 
to each other. We showed the multivaluedness schemat- 
ically by plotting a vs. h/J in Fig. E3 with jagged lines. 

In the intermediate range of (1%), the function a(J/h) 
crosses over from one type of the limiting behavior to 
the other. The numerical data does not seem to suggest 
that (ii) has a universal scaling form of a function of 
a v /(J/h) for all a, J/h < 1. 




FIG. 10: (color online). The dependence of a on h/J as 
given by the condition (7i) = const for different (ii) in the 
chain with L = 300. The lines in the lower panel listed from 
down upward (thin solid, dotted, dashed, long-dashed, and 
dot-dashed) correspond to (I\) = 95, 50, 10, 2, and 1.2, respec- 
tively. The lines in the upper panel listed from down upwards 
(long-dashed, dotted, dashed, and dot-dashed) correspond to 
comparatively small = 2, 1.6, 1.4, and 1.2, respectively. 
The bold lines a = J 2 /IQh 2 and a = J/h display the asymp- 
totic behavior of a for large and small (ii); the line a — J/2h 
corresponds to a = a t h- 



APPENDIX C: NARROW RESONANCES OF 
THE MANY-PARTICLE IPR AS FUNCTION OF 
THE PARAMETER a 

In this Appendix we discuss the positions and widths 
of the narrow peaks of the IPR seen in Fig. 0] 

As we increase a starting from a = 0, pronounced 
peaks of (Je) appear for the difference in the combination 
two-particle on-site energies Se = {e^ + £k 2 ~ £ k 3 — £fc 4 
(jT5|l equal to 

Se pa sah ps J A 

with s = 1,2. They are due to resonant hybridization of 
pairs on neighboring sites (n, n+1) with dissociated pairs 
located on sites (n, n + 3) for s = 1, and (n — 1, n + 2) 
with even n for s = 2, for example. Such hybridization 
corresponds to two single-particle steps by one site, i.e. 
x = 2. 

A specific example for the studied chain with A = 1, 
h/J = 20, and a — 0.05 is the resonance between 
the on-site states |$(3, 4, 6, 7, 8, 9)), |$(4, 5, 6, 7, 8, 9)), 
|^(3,4, 6, 7,8,11)) and |$(4, 5, 6, 7, 8, 11)) (we remind 
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that the arguments of $ indicate the positions of the 
excitations; we have six excitations, and the available 
sites are 1,2,. ..,12). All these states can be obtained 
from each other by moving one excitation by two posi- 
tions. For example, in the first pair the excitation goes 
from site 3, where it has one nearest neighbors, to site 5, 
where it has two neighbors. 

The width of the above peaks 8a can be estimated from 
the condition that the frequency detuning \sah ± JA| is 
of order of the effective hopping integral J A Vk 1 k 2 k 3 k4- 
For s = 1 [an (n,n + 1) «-> (n,n + 3)-type transition] 
the hopping integral is ~ J 3 A/ah 2 from Eqs. (|llfl. (JT^J. 
This gives the width 

6a ~ {J/hf. 

The positions of the peaks a ~ J A / sh and their widths 
are in agreement with the data in both upper and lower 
main panels of Fig. 0] 

For larger a, narrow resonances with respect to a occur 
when 

sa m h w MJA 

with integer s,m,M, and m > 2, M = 1,2. They 
may happen, for example, between pairs (to, n + 1) and 
(n — l,n + 2) with odd n such that n ^ 3k — 1, in 
which case m = 2 and x = 2. A specific exam- 
ple for our chain is the resonance between the on-site 
states |*(1,3,4,6,9, 11)) and [$(1,2,5,6,9,11)) for A = 
1,/i/J = 20, and a = 0.246. Here the excitations on sites 
(3, 4) move to sites (2, 5), and a 2 h ~ J A (in fact, higher- 
order terms in a are essential for fine-tuning the states 
into resonance). In other cases resonances with m > 2 
require more intermediate virtual steps, with x > 4. 

The m > 2-resonances are extremely narrow for ath -C 
a <C 1. For example, for m = 2 their widths are 

Sa < {J/hf/ 2 A 1 / 2 for x = 2, 
fa < J 3 //i 3 A for x = 4. 

In these estimates we used that, from Eqs. (@J, l|12|l 
\V n -in+2nn+i\ ^ J 2 /h 2 for the x = 2-transition (n, n + 
1) <-» (n— 1, n + 2). For the m — 2 and xr = 4-transitions, 
on the other hand, |I4 1 fe 2 fe 3 fe 4 | < J 4 /a 3 h 4 [for exam- 
ple, this estimate applies to a transition (to, n + 1) «-» 



(to, to + 5)]. We note that, from the condition ath -C a 
and the resonance condition sa 2 h — MJA, it follows 
that A 3> J//i, which guarantees the smallness of the 
peak widths. 

Each high-order resonance gives rise to a narrow band 
of resonant a values. All of them refer to a resonant 
transition between the same sites. However, the energy 
difference of these sites is slightly different depending on 
the occupation of remote sites, for example, next nearest 
neighbors. In this latter case, from Eqs. Qllfl. i|12|) . the 
corresponding shift of a is oc ( J//z) 5 / 2 A 1 / 2 . 

A specific example for the studied chain is provided 
by the resonances between two pairs of on-site states, 
[$(2,3,6,7,8,12)) and |$(2, 3, 6, 8, 11, 12)), on the one 
hand, and |$(5, 6, 7, 8, 9, 12)) and |$(5, 6, 8, 9, 11, 12)), on 
the other hand. In both cases the resonant transition is 
fcrmion hopping from site 7 to site 11. Both resonances 
occur for A = l,h/J — 20, but the first corresponds 
to a « 0.2778, whereas the second corresponds to a w 
0.2782. The difference in a comes primarily from the 
different occupation of the next nearest neighbors of sites 
7 and 11. 

The most pronounced peaks in Fig. 0] correspond to 
comparatively small x < 4. However, there are reso- 
nances for higher x as well. The resonating energies have 
to be extremely close to each other for such states to be 
hybridized. The corresponding peaks are very narrow, 
and very high precision is needed to find them numeri- 
cally (sometimes the hybridization appears to be an ar- 
tifact of not sufficiently precise calculations). 

As mentioned above, in the region 0.2 < a < 0.4 the 
positions of the IPR peaks are determined not only by 
the leading-order terms in a, but also by higher-order 
terms. Therefore there are several resonant bands for 
each s,m,M as given by the condition sa m h w MJA. 
This explains why there are several bands in Fig. 0] 
These bands are well separated for sufficiently large h/ J 
and not too large A. On the other hand, for h/J = 10 
and A = 3 the bands of resonances are broadened and 
overlap with each other. When a is not very small 
there emerge also narrow resonances where Se <C J and 
5e/h — * 0. They are responsible for some of the peaks in 
the insets of Fig. 
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